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' This article uses the conformal Einstein equations and the conformal representation of 

spatial infinity introduced by Friedrich to analyse the behaviour of the gravitational field 
\l ■ near null and spatial infinity for the development of initial data which are, in principle, non- 

' conformally flat and time asymmetric. This article is the continuation of the investigation 

^ , started in Class. Quantum Grav. 21 (2004) 5457-5492, where only conformally flat initial 

O^' data sets were considered. For the purposes of this investigation, the conformal metric of the 

initial hypersurface is assumed to have a very particular type of non-smoothness at infinity 
l^i^, in order to allow for the presence of non-Schwarzschildean initial data sets in the class under 

study. The calculation of asymptotic expansions of the development of these initial data sets 
reveals — as in the conformally flat case — the existence of a hierarchy of obstructions to the 
. smoothness of null infinity which are expressible in terms of the initial data. This allows for 

}_( ' the possibility of having spacetimes where future and past null infinity have different degrees 

of smoothness. A conjecture regarding the general structure of the hierarchy of obstructions 
is presented. 

Pacs: 04.20.Ha, 04.20.Ex 



1 Introduction 

This article is the second part of an investigation of the gravitational field near null and spatial 
infinity arising as the time development of initial data sets which are not time symmetric — that 
is, the second fundamental form of the initial data set is non-vanishing. Part I of this study was 
concerned with the analysis of initial data sets for which the 3-metric of the initial hypersurface 
— its first fundamental form — was assumed to be conformally flat |17| . The main result of 
part I was the uncovering of a time asymmetric hierarchy of obstructions to the smoothness 
of null infinity, which if vanishing, would imply that the initial data set under consideration is 
asymptotically Schwarzschildcan to a certain order. The obstructions found in part I are the 
direct generalisations of those found in jJS] for time symmetric, conformally flat initial data. 
In hindsight, the analysis performed in could be readily generalised in two directions: (i) 
removing the conformal flatness assumption, but keeping the time symmetry; (ii) keeping the 
conformal flatness of the initial 3-metric, but taking away the time symmetry; (iii) doing without 
the conformal flatness and time symmetry. The generalisation (i) was carried out in [Tl) , where a 
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hierarchy of obstructions implying asymptotic staticity of the initial data near infinity was found. 
As mentioned earlier, generalisation (ii) was carried out in part I. 

This article is essentially concerned with generalisation (iii). That is, we want to deduce 
asymptotic expansions of the gravitational field near null and spatial infinity for the developments 
of time asymmetric initial data sets which are not conformally flat. From the results in [151 114| 
and part I, we expect that the stationary spacetimes will play a crucial role in our analysis. As 
we are constructing our asymptotic expansions out of Cauchy initial data sets, we are in the need 
of considering a class of initial data big enough to contain strictly stationary data — i.e. not only 
disguised static data. From the results concerning the non-existence of conformally flat initial 
slices in the Kerr and other stationary solutions, it is clear that with conformally flat data, one 
would only be dealing with — more or less sophisticated — Schwarzschildean data ^Hl- Thus, one 
has to move away from conformal flatness. But this is not enough, for in ^ it was shown that 
the conformal metric of Kerrian initial data is actually non-smooth — a feature that is bound to 
pervade all stationary solutions with non- vanishing angular momentum. 

The general strategy of this article can be summarised as follows: given a sufficiently wide 
class of initial data sets, we calculate asymptotic expansions for the initial data of the conformal 
Einstein equations. These expansions are, in turn, fed into the conformal evolution equations 
and used to obtain in a recursive way asymptotic expansions of the development of the initial 
data sets. These expansions describe the behaviour of the resulting spacetime in the region close 
to null and spatial infinity. The recursive nature of the calculations is quite amenable for the 
implementation on a computer algebra (CA) system. Indeed, the crucial results given in this 
article have been calculated using scripts written in the system Maple V. 

As in part I, the asymptotic expansions show generically a certain type of logarithmic non- 
smoothness at the sets where null infinity and spatial infinity "touch" . Because of the hyper- 
bolic nature of the evolution part of the conformal field equations it is quite likely that this 
non-smoothness will be propagated along the generators of null infinity. Hence, the conformal 
completion will be non-smooth. A remarkable feature of our set up is that it allows us to identify 
the parts of the initial data responsible for these logarithmic divergences. As a consequence, 
one obtains a hierarchy of quantities expressible in terms of the initial data — henceforth to be 
referred to as an obstruction — whose vanishing would eliminate a certain type of logarithmic 
divergences. The hierarchy exhibits a subtle structure. Actually, the regularity observed in our 
calculations permits to infer what the general behaviour of the obstructions to any order should 
be. 

In part I, the relative simplicity of the conformal data permitted to show that the vanish- 
ing of the obstructions to a certain order leads to an initial data set which is asymptotically 
Schwarzschildean — but not necessarily time symmetric. In the case of the current article, the 
increased complexity of the initial data sets considered, and of the resulting expansions preclude 
us from deducing an analogous condition with regards to stationarity. That is, if one were able 
to carry the expansions to high enough order one would expect to be able to deduce asymptotic 
stationarity from the hierarchy of obstructions. An analysis of how stationary data and their 
asymptotic expansions fit in the framework of the cylinder at spatial infinity will be pursued 
elsewhere. 

As a continuation to 17 , this article is completely consistent with the notation, conventions 
and nomenclature of part I. The main features of our set up, including the construction of the 
cylinder at infinity and the use of 2-spinors have been discussed there and the reader will be 
duly referred to it whenever it is necessary. The article is structured as follows: in section 2, 
we discuss the existence of solutions to the constraint equations which are expandable in the 
asymptotic region in terms of powers of l/r. The discussion of this section is based in partially 
unpublished results due to S.Dain. As we are using systematically a 2-spinor formalism, in section 
3, we discuss the effects of our class of initial data on the solutions of the structure equations. In 
section 4, it is shown how to fix the conformal gauge in the initial data. This is done by resorting 
to a certain gauge based on conformal geodesies on the initial hypersurface (the cn-gauge). This 
procedure of fixing the conformal gauge is important if one wants to identify the parts of the freely 
specifiable data which are pure gauge. Section 5 is concerned with the expansions of solutions 
of the non-conformally flat momentum constraint. Section 6 does likewise but now with the 
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solutions of the Hamiltonian constraint. Section 7 discusses how a certain regularity condition 

first found by H. Fricdrich can be implemented in our set up. Section 8 is concerned with the 
solutions of the transport equations implied by the conformal Einstein equations on the cylinder 
at spatial infinity from which our asymptotic expansions are deduced. Section 9 speculates about 
the general behaviour of the hierarchy of obstructions and the associated logarithmic divergences 
of the solutions to the transport equations at the sets where null infinity "touches" spatial infinity. 
Finally, section 10 provides some conclusions. 

In order to carry out the calculations here described, a number of assumptions regarding the 
initial data set have been made. These are clearly numbered from 1 to 5. 



1.1 General conventions 

Throughout, we shall fully follow the conventions used in part I. In particular, let {M,gij,v) 

denote a spaeetime arising as the development of some Cauchy initial data (5, haf), Xa[))- Tilded 
quantities will refer to quantities in the physical spaeetime, whereas untilded ones will denote 
generically quantities on an unphysical — i.e. conformally rescaled — spaeetime. The indices n, v, 
A, . . . (second half of the Greek alphabet) are spaeetime indices taking the values 0, . . . , 3; while 
a, /3, 7, . . . are spatial ones with range 1, . . . , 3. The latin indices a, 6, c, . . . will be used in spatial 
expressions which are valid for a particular coordinate system (usually a Cartesian normal one) 
and take the values 1, . . . , 3. The indices, i, j, fc, . . . are spatial frame indices ranging 1, . . . , 3, 
while A, B, C, . . . will be spinorial indices taking the values 0, 1. Because of the use^of spinors, 
the signature of 5^1/ will be taken to be (+,—,—,—), and the 3-dimensional metric /i„/3 will be 
negative definite. 

We shall be concentrating in a particular asymptotically flat end of the manifold S. With- 
out loss of generality, we shall assume that there is only one asymptotic end. Let {y"} denote 
coordinates adapted to that end in the sense that hab = —Sab + C'(l/|?/|) as \y\ 00. We shall 
also assume (5, haf^) to be asymptotically Euclidean and regular in the sense described in part 
I. Accordingly, let {S, ha/3) denote the 3-dimensional, orientable, smooth compact Riemannian 
manifold and i G 5 the point whose suitable punctured neighbourhoods correspond to the asymp- 
totically flat end of {S, hajj)- The manifold S is identified in a standard way with iS \ {i}. Most of 
our discussion shall be carried out in a sufficiently small neighbourhood Ba [i) of radius a centred 
on i. The neighbourhood Ba{i) will be assumed to be geodesieally convex. Unless otherwise 
stated, {a;"} denote some normal coordinates with origin at i based on a /i-orthonormal frame 

A crucial feature of our analysis is that the region of spaeetime near null and spatial infinity 
can be described by means of a certain manifold Ma,K given by 



Ma,K = { (t, q) 



q&ta,p, l-^<T<——\, (1) 

K{q) K{q) J 

where Ca.K = P^^'^Ca and Ca is the blow up of the neighbourhood Ba{i)- The conformal factor 
associated with this representation of the "unphysical" spaeetime is given by 

e = .-^n(i-r^^y (2) 

where is the conformal factor of the initial hypersurface and 

uj = , (3) 

The function k containing the remaining freedom in our setting will be set equal to 

K = p, (4) 

in order to ease our calculations. Accordingly, the critical sets, X^, where null infinity "touches" 
spatial infinity are given by the point where p = and r = ±1. 
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In order to ease the presentation and understanding of the results described in this article, we 
shall be making the following assumption 

Assumption (axial symmetry). The initial data sets {S ^hapjXap) to he considered will he 
assumed to be axially symmetric. 

It follows from well known results — see e.g. ^2] — that the development of the initial data will 
inherit the axial symmetry. It must be stressed that there is no fundamental restriction implied by 
the latter assumption. All the results presented in this article can be extended naturally — aside 
from the associated increase of computational complexity — to the non-axially symmetric case. 



2 Existence of solutions to the constraint equations with 
non-smooth conformal metrics 

As discussed in part I, for technical reasons, we are interested in solutions (/iQ/3,Xa/3) to the 
Einstein constraint equations 

r-X^+Xapx"^ ^0, (5a) 

5"xa/3 - D/sx = 0, (5b) 

which on an asymptotic end of the initial hypersurface S admit expansions of the form 



T-af3 



where h^^^ and xl^p are smooth functions on the sphere. 

We shall perform our analysis not in the physical spacetime which arises as the time develop- 
ment of (5, /iq,/3, Xq/3)j but in a conformally rescaled version thereof — the unphysical spacetime. 
Writing as usual 

= n^g^,u, (7) 
we make the following assumption (which was also made in part I): 

Assumption 1 (maximal initial data). It shall he assumed that 

S = 0, Vlx^x^O onS, (8) 
where E denotes the derivative of along the future g^unit normal of S. 



Under the latter assumption, the constraint equations (|5a|l - (j5b|l imply on the unphysical 
manifold S the equations 

[d^D" - irj z9 = with d = (9a) 

D'^ {i^^p) = 0, (9b) 
where D is the connection associated with the conformal 3-metric 

hap = ^'^hajS, (10) 

and the 3-tensor tpap is given by 

V'a/3 = f^"^Xa/3 = 1?^Xq/3- (H) 

The equations (|9a|l and Ij9bp are to be supplemented with the asymptotic (boundary) conditions 

\x\d^l, ^p^p^of-^), as|x|^0, (12) 
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where {a;"} are normal coordinates with origin at i. 

Conditions for the existence of solutions to the equations 19a|l - (|9b(l with the above boundary 
conditions, and such that the corresponding physical fields have the asymptotic behaviour © 
have been given in 171 under the assumption that the conformal metric haf^ is smooth in a 
neighbourhood Ba{i) of infinity with respect to some normal coordinates. Unfortunately, it seems 
that the assumption of smoothness of the conformal metric is not enough to discuss stationary 
spacetimes. In ^ it has been shown that for Kerrian data derived from the Kerr spacetime in 
Boyer-Lindquist coordinates, 

hab^hii^ + \x\'hii'\ (13) 

where /ij^P and h'^j^^ are analytic in Ba{i)- The important fact to be noted here is that \x\^ € 
C'^'"{Ba{i)), whence also hab & C^'°'{Ba{i})- In the light of the results given in [31^] it is very 
unlikely that there are slices in the Kerr solution with a better smoothness at infinity. Now, 
stationary spacetimes are, to leading order, Kerrian 1 , and thus, similar non-smoothness is to 
be expected from any other stationary solution with non-vanishing angular momentum. 

Existence results for solutions to the equations (|9a|l and Ij9b|l under assumption of a conformal 
metric of the form given in equation (|13() and which are expandable in powers of |x| are only avail- 
able under the premise of axial symmetry |^|2j. Nevertheless, results for the non-axisymmetric 
case seem plausible. 

Dain's existence results use the fact that for axially symmetric initial data sets it is possible 
to explicitly write solutions to the momentum constraint. These solutions are calculated from 
freely a specifiable function which acts as a potential. So, if the metric hab has an axial Killing 
vector ?7", then a solution of the momentum constraint, equation l)9b|l is given by 

where 

r = , ^, e'^'^r^bDciu, Cr,Lo = 0. (15) 

hdeV V 

The potential u in the last equation is calculated from a spin real function A and a constant 
via 

uj = g^Asin^ + J^(-3 cos 61 + cos^ 9), (16) 

where 3 is the standard NP "eth" operator. 

It turns out that solutions to the momentum constraint constructed by the aforediscussed 
procedure satisfy the regularity condition given by Dain & Friedrich in [Jj. Indeed, 

Theorem 1 (Dain, 2003). Assume that a metric of the form ^13\) has an axial Killing vector 
r]". Let ^p"'' be given by \T^, and |2B). //|a;|A € E°"{Ba{i)) then \x\^i}abil'''^ € E°°{Baii)). 

The purpose of the condition \x\^ipabip°''' S E°°{Ba{i)) is to ensure that the solutions to the 
Hamiltonian constraints are free of In \x\ terms. This is the content of the following 

Theorem 2. Let hab be a metric that on Ba{i) £ S is of the form 1^1S\) with x'^h^^j^^ = 0, which 
is furthermore smooth on S — S \ {i}. Let ipab be given by equations fl5\) . and 1^16]) . Then 

there exists a solution i) to fffal) which is positive, satisfies the boundary conditions and on 

Ba{i) has the form 

^eE^{Ba{i)), ^9(z) = l. (17) 

\x\ 

3 The structure equations on S 

Let = v^Cq, where Bq is the future g-unit normal of 5, and let SU{S) be the bundle of space- 
spinors associated with , the spinorial counterpart of . In part I, a certain submanifold 
Ca C SU{S) given by 

Ca = |<5(p,t) e SU{S) I IpI < a, t e 5C/(2,C)1, (18) 
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was defined. The manifold Ca inherits from SU{S) a 5u(2, C)-valued connection form Cj^^ com- 
patible with the metric hap and 1-form , the solder form, such that 



hapdx'^ ® dx^ = h ABC DC 



where 



AB 



.AS, 



Let CCD be the vector fields dual to the forms , in the sense that 

(ct^^,ccd) = h^^cD^ CCD = Ccjjdp + c^j^X+ +c^^X^, 

where 

1 + 1 + _ 1 _ 

Cab = XAB, C^B = -ZAB + C^B, C^B = + C^B" 

General considerations similar to those in 8 require that, 



(19) 
(20) 

(21) 
(22) 



C^B = Y^VAB + Z^ZAB, 

The relation (|19|l shows that a knowledge of the conformal metric hap is equivalent to a 
knowledge of the solder form , and this, in turn, is equivalent to a knowledge of the frame 
vector field cab- Hence, and because of the systematic use of space spinors in this work, we 
shall regard the regular parts, c^bj of the frame coefficients, c^g, as encoding the non-conformal 
flatness of the conformal metric hap- 

Now, assuming that the conformal 3- metric, hab is of the form (|13|) . it follows from the relations 
(|19ll - (|22ll and a "bookkeeping" argument that 



C-AB — \C )ab + P\C JAB ' 



(23) 



where {c^)ab ~ ^(P^) ^ ^(p) analytical spinorial fields on Ca- More precisely, 

in order to have a class of initial data sets large enough to contain stationary data we make the 
following assumption: 

Assumption 2. The frame fields cab on Ca are such that 



-AB 



in'Ak+p(.n 



AB 



where {c^)ab ^ ^(P^) '^'^^ ic^)AB ^ ^ip) ^ire analytical spinorial fields on Ca- 

Remarks. The spin weights of the diverse components of imply under assumption |21 the 
following expansion Ansatze: 



oc p 2q CO p 2q 

^"^ = E E E >;^2,,™T,/v,p^ - E E E p^ (24a) 

p=2 q=2 m=0 p=2 q=0 m=0 

oc p 2q oc p 2q 

= E E E Z;,,q,nT,^"'q.2P', = E E E Z,.,2q,rnT2q"'qP', (24b) 

p=2 q=2 m=0 p=2 q=0 m=0 

where the coefficients Y^2q m ^'^'^ ^p-2q m '"^^^ complex numbers satisfying the reality conditions 



Z. 



p\2q,m \ J ^ p;2g,2g— m 



V~ — \m+q 7+ 

1 ^p:2q,m — \ -'-j '^p: 



2q,2q—m - 



(25) 



Note that, for example 



Y- 



2! 



^2;0,0-'0 I ^2;4,m^4 
2 

^2;2,m-'^2 1 



^ +3! 



m=0 



E ^;2,m^2'"l + E ^3;6,m^f 
m—O m 

4 

^3;0,0^0 + E/ ^3,4,m^' 



6 3 



m=0 
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so that the analytical function contained within the second pair of round brackets is responsible 
for the non-analyticity of the frame vectors — and of that of the conformal metric. Similar 
arguments can be made with the remaining coefficients Y~^, and . More precisely, one has 
the following 

Lemma 1. // the frame fields satisfy assumption\^ then in Ba(i) the conformal metric is of the 
form 



ha, = h[[^ + \x\hii'\ 



with h^^^ and h^^^^^ analytic. 

Thus, assumption [3 implies a conformal metric which is rougher at infinity than that of the 
Kerr data discussed by S. Dain in jj- 

Finally, it is noted that under the assumption of axial symmetry (assumption 0) the Ansatze 
(I24a|l - (j24b|) simphfy to 

oo p oo p 

^+-EE^pk,^2,V2P^ y-^T.T.Ypi^.Mp", (26a) 

oc p oc p 

^+ = E E K;2,,,T,^V2P'', = E E ^pIs,,,^^, v^ (26b) 

p=2 9=2 p=2 q=0 



while the reality condition H25() becomes 

^p;2g,g ^ p;2g,gi ^p-,2q,q ^p;2q,q- 



^V:2q.q ~ y p\2q.qJ ^p^q.q ~ ^p;2q,q- (27) 



3.1 Expansions for the connection coefficients 

Now, recall that the contraction of the connection form Lu"^g with the frame cab defines the 
connection coefficients ^cdab- More precisely, we write 

ICDAB = {i^^BiCcd)£EA = Tri^ACXBD + (-BDXAc) + ICDAB = -IcDAB + ICDAB, (28) 

2p p 

where 

loiCD = 0. (29) 

Moreover, it can be checked that 

71100 — ~7ooiii 70010 = 70001: 7iiio = 7iioi- (30) 

As seen in part I, the regular part of the connection coefficients, Jcd^b vanishes if the confor- 
mal metric is fiat. However, if the conformal metric satsifies assumption |21 then the connection 
coefficients depend on the frame coefficients c^g . In order to calculate the connection coefficients 
we make use of the commutator equations 

E ^ F „, s . -F^„„„ , 



[cAi3,Cci3] = 2|^7^S (C^D) ~^CD (A^B) j ^EF + [Cj^bCcD ~ ^CD^^'Ab) ^- (31) 

The latter imply the radial equations, 

-^9p{pc\^) = "YAAooipcti) - 7AAii(pcJo) - -i=7AAii, (32a) 

-^dp{pc^^) = jAAOoipCii) - lAAllipCoo) - -^JAAOO- (32b) 

The latter equations can be used to obtain expansions for the connection coefficients 70000 j 
70011 = —71100 and 71111 in terms of the coefficients in the expansions of y+, Y~ , and Z~ . They 
also show that the frame component Z^ is not independent, but can also be expressed in terms 
of Y~^, Y~ , and Z~ . Furthermore, one has the following 
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Lemma 2. Under assumption\^ it follows that the connection coefficents 7000O; 7ooii; 7iioo o^i^-d 
71111 such that 

lABCD = {jabcd)^^^ + p{iabcd)'^^^\ (33) 

where {'^abcd)'^^^ = 0{p^) and {^abcd)'^^^^ = ^{p) '"^g analytic fields on Ca- Furthermore, the 
reality condition {'^abcd)^ ~ —Iabcd is automatically satisfied if the conditions \25l hold. 

The proof of this lemma foUows again from a "bookkeeping" argument. Note that the con- 
nection coefficients 70001 and 71110 are not fixed by the procedure we have discussed above. We 
shall come back to them later. 



3.2 The curvature spinors 

In order to determine expansions for diverse components of the curvature we shall make use of 
the structure equations. Evaluating the second structure equation — equation (36b) in part I — 
on CCD cef one obtains the following radial equations 

^ a - ^ 1 /- - ^ 1 - 

—^OplOOAB + - i lOOOOZAB ~ 'JOOllVAB H ^700AB 

= 70000711AB - 70011700AS - -^SABoo - "^^^y^^' (34a) 

In- If- - 1 - 

— ^Cp7iiAS + - \ ^IWQZab ~ lllllVAB H — ;=7llAB 

= 71100711AB - 71111700AS - -j^SABii - ^yl^^^-s, (34b) 

where sabcd denotes the spinorial representation of the tracefree part of the Ricci tensor of hap, 
that is, 



sabcd = crliB'^cD^i^j yafi - -^rhap j , (35) 
where r^p is the Ricci tensor of hap and r — h°'^rap- Because of its symmetries, one can write 

SABCD = SoeABCD + ^l^^ABCD + ^^^ABCD + ■'^3'^ABCD + ^i'^ABCD- (36) 



The equations (|34all and (|34b|) have been used to obtain expansions for the components of sabcd 
in terms of F+, , and . More precisely, one has the following 

Lemma 3. Under assumption\^ one has that the components of the Ricci spinor are of the form 

s," = s) ' + -s) ' 

3 p 3 

where s^P ~ 0{p) a nd s^ "* — 0{p^) are analytic fields onCa- The reality condition {sabcd)'^ — 
SABCD is satisfied if \2fA holds. 



The tracefree part of the Ricci tensor and the Ricci curvature scalar are not independent from 
each other, but related via the Bianchi identity 

D"^^ SABCD = Idcdt. 

6 

The latter identity implies three equations which can be used to write the Ricci scalar r and the 
connection coefficients 70001 and 71110 in terms of the frame freely specifiable data — i.e. , 
Y~ , and Z^. It follows that 
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Lemma 4. Under assumption 2 it follows that the Ricci scalar is of the form 

P 

where r^^^ — 0{p) and r*^^^^ = 0{p^) are analytic fields on Ca- Furthermore, the components 
70001 cw^c^ 71110 of the connection are such that 

^ABCD = {'^ABCd)^^^ + p{iABCd)'^^^\ 

where {'^abcdY^^ — ^{p^) {'^abcd)'^^^^ = ^{p) ^.'"'^ analytic fields on Ca- Furthermore, the 
reality condition {'^abcd)^ = ^Iabcd is automatically satisfied if the conditions \25\ hold. 



4 Fixing the conformal gauge on S 

In our set up, the functions F+, , and determining the frame — via the relations 1)21(1 
and are used to specify the part of the free data which is more commonly encoded in the 

metric. In their expansions, these functions contain terms which are pure gauge. The reason for 
this is the following: given ■& and ipa/S solutions of the constraint equations H9al) and (|9b|l and a 
positive function, on S, it is well known that the transition 

yields another solution to the Einstein constraint equations with the same physical data. In order 
to fix this conformal freedom we shall use a certain gauge based on spacelike conformal geodesies 
starting at i. We shall refer to this gauge as the cn- gauge. 

The cn-gauge is defined as follows: consider the 3-dimensional conformal geodesic equations, 

xf^Wpi" = -2b0x'^x°' + xpx'^b", (37) 

X^Vpba = bpX^ba - ^bpb'^Xa + (^Sal3 + -^rhafi^ x^ , (38) 

where a:"(i) is a curve on 5, and 5" an associated 3-dimensional 1-form. We supplement the 
conformal equations with the initial conditions, 

x{Q)^i, Xf}x^ = -l, 6(0) = 0. (39) 

If a is chosen small enough, there exists a unique solution to these equations on Ba{i). Further- 
more, there exists in Ba{i) an unique conformal rescaling such that, 

bpx^ = on Ba{i), (40) 

can always be found. 

Definition. A metric in the conformal class for which the condition \4-U\j is satisfied along the 
solutions of the 3-dimensional conformal geodesic equations \31\) and will be said to be in 
the cn-gauge. 

It is noticed that for conformally flat data, if x is the tangent to (standard) geodesies starting 
at i with XjjX^ = — 1, and if one requires 6 = in Ba{i), one is automatically in the cn-gauge. 

Let Xab and Bab be, respectively, the space spinors corresponding to the 3-dimensional 
vector and the 1-form These can be written as, 

Xab = XxXAB + XyyAB + X^zab, Bab = B^xab + ByyAB + B^zab- (41) 

Using these two spinors and the spatial Infeld symbols, one can produce a spinorial version of the 
conformal geodesic equations H37|) and (|38|l . 
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Following our overall strategy, one can attempt tosolve the 3-diniensional conformal geodesic 
equations together with the cn-gauge assumption by means of formal expansions. Spin weight 
considerations plus the initial data for the conformal geodesies show that 

X, = l + 0{p), Xy^O{p), X,=0{p). (42) 

Similarly, for the components of the spinor Bab one has that 

B.,^0{p), By^O{p), B, = 0{p). (43) 

Once the conformal geodesic equations have been solved, requiring the cn-gauge condition x°'ba — 
e^'^' e^^ XcdBab — renders further relations between the coefficients in the expansions of the 
frame components. Indeed, an explicit calculation with Maple V shows that 

Lemma 5. Under assumption\^the cn-gauges imply that the coefficients Y~2q m with 2 < p < 7, 
^^(l^P,^^Tn<2q are determined by the algebraic combinations of the coefficients Y^2q m 
and Z~.2g ,„ with 2<p<7,0<q<p,0<m<2qof and . The only exceptions are the 
coefficients Y^2 m? m = 0, 1, 2 and Y^q g which remain freely specifiable. 

The full details, being not very illuminating, shall not be given here. Further consequences of 
the cn-gauge are that in the decompositions 

^ = + s,^s^'^ + -s^"\ (44) 

p ■' p 

with j = 0, . . . , 4 the "smooth parts" are of the form 

rW=0(p2), =0(p2), (45) 

instead of being 0{p) — cfr. with lemma 4.7 in 8 . The "non-smooth parts" do not change order. 
It may also be of interest to note that if the cn-gauge condition holds then 

X, = 0(p3), Xy^O{p^), X,^0{p'), (46) 

and 

B,^0{p'), By^Oip'), B,^0{p^). (47) 



5 Expansions of the momentum constraint 

We are now faced with the problem of choosing a class of solutions of the momentum constraint, 
equation (joEJ, which is large enough to encompass the second fundamental forms of stationary 
solutions with non-vanishing ADM angular momentum. In agreement with the existence results 
given in section |2 we shall further require that the solutions satisfy the regularity condition 
ipapip"^ & E°°{Ba{i))- The latter requirement is essentially satisfied if one considers solutions to 
the momentum constraint without linear momentum. Our selected class of solutions will be, in 
turn, lifted to the manifold Co, where formal expansions in terms of powers of p will be calculated. 

The results by S. Dain discussed in section 2 strongly use the fact that the initial data pos- 
sesses axial symmetry to explicitly calculate — via a "potential" , A — solutions to the momentum 
constraint Ij9b|l . However, with future applications in mind, we are interested in a procedure that 
can be easily generalised to a non-axially symmetric setting. 

For the time being let us write 

■00/3 = -00/3 + V'a/?, (48) 

where both -0^^ and tpa/i are symmetric, tracefree tensors. The former is a "seed" tensor and 
the latter an unknown to be determined upon substitution in the momentum constraint. The 
standard strategy at this point is to write ijjap in terms of the Killing operator of a 1-form, Va 



10 



— the York splitting. For reasons to be explained in section 7, we will not proceed in this way, 
but will solve directly for ipafs — see, however, the remarks in subsection 15 .11 

In order to make contact with the conformally flat case studied in part I, we shall consider as 
our "seed" tracefree tensor ip'^^^ the solution to the conformally flat momentum constraint used 
in part I — see assumption 3 in jl7|. As in part I, let 

^2 = (49a) 
= ^t^ = = 0, (49b) 

^? = ^{Q2+ ^Ql)T^^ - ^QsT^'o ~^iQ2~ iQi)T^\, (50a) 
= ^(^2 + iQi)TA - ^QsT^'i ~ ^iQ2 - ^Qi)TA (50b) 

= "4(^2 - ^Qi)TA - ^QsTA + 4(Q2 + *Qi)r2°2, (50c) 
= - 0, (50d) 



tl = ^i-Ji + ^J2)T^% + ^^J^T^'o ^J2)T^%^ (51a) 

p6 pi p6 

= (51b) 
^3 = ^{Ji + ^J2)TA ~ ^^^3^/2 - ^(-Ji + iJ2)T^\, (51c) 
= = 0, (51d) 

were A, Ji, J2, •^37 Qi, Q2, Qz G Recall that assumption — axial symmetry — requires 
Q\ = Q2 = Ji = J2 0. Furthermore, we write 

= ^ E E i.,n-4;2,,,r2, V^+.p", J = 0, . . . , 4 (52) 

^ n=2g=2 

where 

((i4,n-4;2g,9 - io,«-4;2<7,q)(g + l)g + 4Lo,„_4;2g^g(n + 3)2) r- — — — ■ 

= (n + 3)(2(„ + 4)(„ + 2)-(<, + 2)(<,-l)) ^(9 + 2)(« - D. (53a) 

((Lo,n-4;2g,9 - i4,n-4;2g,g)('? + l)g + 4L4,ri-4;29,9("- + 3)2) — — — 

^^--^^^-^ (n + 3)(2(n + 4)(. + 2)-(, + 2)(,-l)) ^(^^^K^, (53c) 

and n = 2, 3,4, . . ., axial symmetry (assumption 0) already being assumed here. The coefhcents 
Lo;2q,q and I'4;2g,g are arbitrary complex numbers satisfying the reality condition 

Lo;2q,q = -^^4;2ij,5- (54) 

We shall make the following 

Assumption 3. The totally symmetric spinor tp'j^g^j-, associated with the "seed" tracefree tensor 
ip'^p in the York splitting is assumed to be of the form 

i'ABCD = i'ABCD + i^ABCD + '^ABCD + ''PaBCD: (55) 

where 

'^ABCD = %'^ABCD + '^I'^ABCD + '^2^\bCD + '^I'^ABCD + V'l^ABCU' (56) 

and i denotes any of the labels A, Q, J , A. 
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The symmetric, tracefree tensor ipap has a spinorial counterpart iPabcd of the form 

■IpABCD = ■IPo'^ABCD + i'lAsCD + ^2<^ABCD + ^I'^ABCD + i'i'^ABCD- (57) 

Under assumption |31 one has that 

D^^'^'abcd = 0{p). (58) 

Consistent with our general strategy of finding solutions to the constraint equations which are 
expandable in powers of p near infinity, we put forward the Ansatz 

oc p 2q 

'^J^Y^ H V'j,p;2g,mT2g" (59) 

p=2 i3=|2-i| "1=0 

where the coefficients ipj,p-2q,m satisfy the required reality conditions — see part I. Under assump- 
tion — axial symmetry — our Ansatz simplifies to 

oo p 

'l'3,P;2q,qT2q%-2+JP^- (60) 

p=2q=\2-l\ 

In this case, substitution into the momentum constraint yields algebraic equations for the coef- 
ficients 'ipj,p;2q,q- These allow us to write the coefficients V'i,p;2q,g, ip2,p;2q,q and ip3,p:2q,q in terms 
of the freely specifiable data and in terms of the coefficients '4'o,p;2q,q and '4'i,p;2q,q- The latter are 
related to each other by the reality condition 

'$4,p;2q,q = '$0,p;2q..q- (61) 

The relations discussed in the above paragraph have been explicitly calculated in the CA system 
Maple V up to order p = 4 in the Ansatz H6U|I . 

5.1 A remark concerning the York splitting 

For the sake of completeness, it is of interest to analyse what happens if instead of using the 
procedure described in the previous subsection, one fully uses the York splitting. Accordingly, let 
us write 

= V'a/3 + (^w)a/3, (62) 

where 

2 

{^v)ap = DaVp + DpVa ~ -ha^D^v'' (63) 

is the conformal Killing operator. Let vab be the spinorial counterpart of the 1-form Va on which 
the conformal Killing operator in the York splitting — equation (|62|l acts, and write 

VAB = V^XAB + VyyAB + VzZAB, (64) 

where the coefficients Vx, Vy and satisfy the reality conditions 

Vx^Vx, Vy^v^. (65) 

Under assumption 13 one has that 

D^^'i^ABCD = 0{p). (66) 
The latter suggests the following Ansatz for the components of vab- 

oo p 2q 

'"x = YYY '"^,p;2q,mT2g"'gpP, (67a) 

p—1 q—0 m— 

oo p 2q 

= YYY ^y,p;29,m?29 " -n/^. (67b) 

p—1 q—1 m—0 

oo p 2q 

^YYY Vz,p;2q,mT2g"'g_ipP, (67c) 
p—1 q—1 m—0 
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where because of the reality conditions one has 

Vy,p-2q,m = (-l)''+"+'lJ.,p;2,,m. (68) 

Under assumption — axial symmetry — the above expressions reduce to 

oo p 

p=l q=0 
oo p 

p=l q=l 
oo p 

=^^V^,P;2q.qT2q"'q-lP^: 
p=l q=l 

with 

^y,p-.2q,m — '^z,p;2q,m- 

The substitution of this Ansatz in the spinorial counterpart of the 
{^v)abcd, and then in turn in 

D'''' [i^ABCD + {^v)abCD^ = 0, 

renders a series of algebraic equations for the coefhcients Vx.p;2q.q, Vy^p-2q,q and Vz^p-2q,q- Calcu- 
lations using Maple V reveal that at least for p — 1, ... ,4, the coefficients Vx,p-,2p,p, '^y,p;2p,p and 
Vz,p;2p,p — i-e. for the sector q = p — can be written in terms of the freely specifiable data — i.e. 
the coefficients Yj^2q,q^ ^p;2q,q^ io,p;2g,g, i4,p;29,g, and the constants A, J, Q3. Now, for q < p 
the situation is different. It is found that for arbitrary specifiable data, there is in general, no 
solution to the algebraic equations linking the coefficients Vx^p-2q,q, Vy^p-2q.q and Vz.p-2q,q- Never- 
theless, conditions for the existence of solutions can be found. These restrictions on the specifiable 
data can be understood by recalling that generic solutions of the momentum constraint, equation 
(|9b|) . will have a much more complicated asymptotic behaviour than the one we have imposed 
through the Ansatze Ht)7a|) - (|67c|) . One could have solutions containing logarithmic terms, for 
example. Consequently, conditions on the free data are needed in order to have consistency with 
the asymptotic behaviour which is being pursued. The conditions are given in the following 

Theorem 3. Under assumptions 1, 2 and 3, necessary conditions for the existence 0/ solutions 
to the momentum constraint, equation iff&j) . which on Ca are expandable in powers of p up to 
order 0{p'^) are 

^p-M.q^^-^y^^^^^v-M.q^ p = 2,...,5, g = 2,...,p, m = 0,...,2g, 

A further peculiarity of solutions to equation H5.1|l under the Ansatze (|67a|) - (|67c() is that even 
if the conditions given in theorem O are satisfied, the coefficients W2:.p;2p,gi 'Vy,p;2p,q and Vz^p-,2p,q 
are not fully determined by the algebraic conditions implied by the substitution of (|67a|l - H67c|l in 
equation 1)5. l|l . More precisely, if the conditions given in theorem O hold then it follows that 

Vy,p;2q,q + {^^Y^"^^^ V z,p;2q,q — Cp-2q,q (71) 

for p = 1, . . . , 4, q = 1, . . . ,p and where Cp-2q,q denotes some expression depending on the freely 
specifiable data. Recall that besides, the reality condition H68|l must also be satisfied. The 
reason for this behaviour can be understood as follows. It is well known that the idea behind 
the York splitting, equation l|()2|) . is to obtain an elliptic equation for the 1-form out of the 
underdetermined system l|9bp . Now, solutions to a system of elliptic equations contain global 
information about the underlying initial manifold S. This information is not captured if one 
attempts to calculate formal solutions in a neighbourhood Ba(i) of infinity. This same feature 
pervades the approach for solving the momentum constraint followed in the previous section. 
Moreover, a similar phenomenon can be observed when "solving" the Hamiltonian constraint by 
the same methods — see in particular the discussion in part I. 
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(69a) 
(69b) 
(69c) 



(70) 

conformal Killing operator 



6 Expansions of the Hamiltonian constraint 



As discussed in part I, the Licnerowicz equation, being a scalar equation is easily translated into 
space spinor language. Indeed, one has 

(^D^^Dab -^r^^ ^ VabcdV^-^^^^^-^. (72) 

In order to construct formal expansions of the solutions of this equations on Ba{i) — respectively 
Ca — we adopt the local parametrisation 

^=^^+W. (73) 

The term C//|x| (the Green function) contains information regarding the local geometry near 
infinity. In the conformally flat case on has U = 1. On the other hand, the term W contains 
information of a global nature. For example, W{i) — m/2, where m is the ADM mass of the 
initial data set (5, hap, Xaf}) — the global quantity par excellence. 



6.1 The construction of the Green function 

In order to calculate the function U, we shall make use of Hadamard's parametrix construction 
— see e.g. jSIEl- The term J7/|a;|, being the Green's function of the Yamabe operator — see 
below — , corresponds to a solution of the equation 

D^D^ - Ir] (^] = inS,, (74) 



where 6j is the Dirac distribution centred in i. Now, consider a function U of the form 



f/ = ^f/pN^ (75) 

p=0 

where Up = Up{\x\), that is, the coefficients in the Ansatz are still allowed to have p-dependence. 
Under this Ansatz, equation (|74|l implies the following hierarchy of ordinary differential equations: 

D-^p'D^Uo = -]^{D'^Dy + 6)C/o, U{i) - 1, 

D'^p^D^Up = -^{D'^D^p^ + 6 - Ap)Up - j^Lh[Up-il p - 1, 2, . . . 



where 



L^,[f] = {D'^D^ - l/8r)/. 



denotes the Yamabe operator of the metric hap applied to a smooth function /. The solutions 
to the above equations can be given recursively by 

/If ds\ 
Uo = exp J {D'^DaS^ + 6)— ) , 



s J 

Uo r Lh[Up] sP 

{4p-2)pP+^Jo Uo 



Up+i = -— _„4-i / 77^ ds, p = 0,l,. 



Now, given a conformal 3-metric hap, the above procedure does not necessarily render a conver- 
gent U as given by the series (I75|l . One needs extra assumptions on h. For example, if h^p is 
analytic in a neighbourhood of i then U given by H75|l is also analytic. 

Under the more general requirements of assumption |21 the best one can hope for is a function 
[/ of a similar form — i.e. U = U^^^ + pU^'^\ The latter will be the case only if further conditions 
on the conformal metric hold, for some experimentation reveals that free data of the form given 
by assumption 121 render coefficients Up with logarithmic terms in \x\. Direct calculations with 
Maple V yield the following 
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Theorem 4. Under assumption\^ the coefficients Uq, . . . have a logarithm- free dependence 
on \x\ if and only if 





- -y+ 

^ -'^2:4,2' 






^3;0,0 


= 0, ^^2,1 


= 0, 




^3;6,3 


^ -^3:6, 3' 






^4;4,2 


- -y+ 

— -'^4;4,2! 


^4;8,4 = 


-y+ 

4:8,4' 


^5;6,3 


- -y+ 

^ -^5:6, 3' 


^5:10,5 


-y+ 

^5;10,5 



// t/ie above conditions hold, then 



Uo 


= i + o(N') 












= 0(1x1°), 




= 0(1x1-1), 


U5 


^Oi\x\-'). 



Remark. The conditions given in theorem 0] allow to calculate an expansion for U up to 
order 0{p^). 



6.2 Construction of expansions of W 

The substitution of the parametrisation (|73|l in the Hamiltonian constraint, equation (|9a|l . renders 
the following equation for the function W in Ca'- 

D^'^Dab ~lr)w^ l^ABCD^^''''''{U/p + W)-\ (78) 



Under the assumption of axial symmetry (assumption 0), we shall be looking for solutions on Ca 
of the form 

oo p 

W^m/2 + Y,Yl ^^^p;2...T^2,>^ (79) 

p=l q=0 ^' 

with 'Wp-2q,m S M SO that is a real function on Ca- 

As in the conformally flat case discussed in part I, the substitution of the Ansatz H79|l into 
equations H78|l allows to determine the coefficients Wp-2q,q, with q < p. The coefficients Wp-2p,p 
— like the ADM mass itself — are, on the other hand, not fixed by this procedure, and shall 
remain free in our construction. Note, however, that in an actual solution of the constraints, 
these coefficients are functions of the free data. This dependence can be analysed, for example, 
by writing an integral representation of the solution — see e.g. the discussion in part I or in [S]. 

In part I, it was shown that the dipolar terms in the expansions of W can be, without loss of 
generality, chosen to be zero. This could be done because our setup still allowed for a translational 
freedom. In our case, due to the use of cn-gauge to fix the conformal gauge this is in general not 
possible. Nevertheless, in order to ease the complexity of our computations we make the following 

Assumption 4. We shall restrict our assumption to initial data sets such that the function W 
in the parametrisation d — U / p + W of the conformal factor is of the form 

w^=^+o(ixn. 

Assumption ^ implies that in the expansion H79I) one has wi;2,m = for m = 0, 1, 2. For the 
purposes of this article the expansions of W have been explicitly calculated up to order O(p^). 
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7 Regularity conditions on the initial data 

As in part I, let us expand the components of the Weyl spinor (j)ABCD as 

oo ^ p 2q 
p=\2-]\ ^' q=\2-j\ m=0 

where the (complex) coefficients aj^p-2q,m have t dependence — that is, aj,p-2q.m — 0'j.p;2q,m{T)- 
Friedrich — see — has shown that the solutions to the transport equations implied by the 
conformal Einstein equations at spatial infinity develop a certain type of logarithmic singularity 
at the sets where null infinity "touches" spatial infinity unless the initial data is such that 

ao,p;2p,fc(0) = a4,p;2p,fe(0) (80) 

with p = 2, 3, . . ., and k — 0, . . . ,2q holds. 

In part I, the regularity condition (|80|l was rephrased in terms of the freely specifiable terms 
available in the case of conformally fiat initial data sets. It essentially implied that the coefficients 
Lo,p-4;2p,fe and -L4,p-4;2p.fc encoding part of the higher multipole structure of the conformally fiat 
initial data sets have to be alternatively real or pure imaginary numbers depending on the parity 
of k — see theorem 3 in jTJj. 

In the case of initial data of the form being considered here, we have the following 

Theorem 5. For an initial data set satisfying the assumptions 1-4, the following conditions are 
equivalent: 

(i) 4'o,p;2p,p{0) = 04,p;2p,p(O), p = 2, 3, 4, 5, 6; 
(n) 



y+ - 

2:4,2 — 


2 (^4,-2:4,2 " 


f ^0,-2;4,2) 1 


(81a) 


Y+ - 

^3:6,3 ~ 


^ (^4,-1:6, 3 ^ 


^" io,-l;6,3) : 


(81b) 


Y+ - 

-"^4:8,4 ~ 


96 

-p- (io,0;8,4 - 

5 


^' i4,0;8,4) , 


(81c) 




: —100 (io,l;10,5 + i4,l;10,5) : 


(81d) 


-'6:12,6 ~ 


1080 , ^ ^ , 

^ (-t^0,2;12,6 + J^A,2:12fi) ■ 


(81e) 



This result naturally leads to our final assumption. 
Assumption 5. The second fundamental form satisfies the condition (ii) in theorem\^ 

8 Solutions to the transport equations and obstructions to 
the smoothness of null infinity 

We shall now proceed, in a parallel way as it was done in part I, to discuss the properties of the 
solutions of the transport equations implied by the conformal Einstein equations — see equations 
(131), (132) and (133) — in the case of initial data sets satisfying assumptions 0-5. For details 
on the structure of the transport equations we remit the reader to part I. We just recall that the 
transport equations implied by the conformal Einstein equations upon evaluation on the cylinder 
at spatial infinity naturally group in two subsystems. One, called the v subsystem, contains the 
evolution equations for the frame, connection, and the components of the Ricci tensor of the 
Weyl connection. The other one, the so-called Bianchi subsystem, provides equations for the 
components of the Weyl tensor. As in references |15II14[[T7| . we shall be restricting our attention 
to the logarithmic terms arising in the components of the Weyl spinor 4>abcd — the reasons for 
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this way of proceeding have been discussed in the aforementioned references. To this end we 
recall that we are decomposing the components of the Weyl tensor as 

oo 

P=|2-i| ^ 

where 

E il'^M.Ar)T,^V,^y (83) 

q=\2-j\ k=0 

Retaking a notation already used in part I, we use Q(t) to denote a generic polynomial in r, 
while VkiT) will denote a generic polynomial of degree fc in r such that 7'fc(±l) = 0. It will be 
understood that the Q(t) and Vk{T) appearing in different equations and/or components are also 
different from each other. 



8.1 Lower order solutions 

Our first result is the following 

Theorem 6. Under assumptions 0-5, the solutions of both the v and Bianchi transport equations 
for p = 0,1,2 are polynomial in r. In addition, the solutions of the v transport equations are 
polynomial for p = 3. 

A first difference with respect to the cases analysed in paper I and also in [151 114L [TT] , is that 
extra conditions are needed to calculate the solutions of the p = 3 Bianchi transport subsystem. 
Namely, one has that 

Theorem 7. Under assumptions 1-5, the Bianchi transport constraint equations are satisfied if 
and only if 

i4,-2;4,2 = — io,-2;4,2. (84) 

Furthermore, if the latter condition holds then the solutions of the Bianchi transport equations 
for p — 3 are polynomial in r. 

Remark. It is noted that the condition H84(l arose in paper I as a regularity condition which 
precluded the appearance of logarithms associated to the "highest possible" harmonics. In the 
present case (|84|) implies, in turn, that 

^2t4,2 ^ ^;4,2 ^ ^2-A,2 ^ ^2:4,2 ^ 

from where it follows that c^g ~ 0{p^). 



8.2 The solutions at order p = A 

Now, let the assumptions 1-5 together with the condition H84|l hold. The CA computations show 
that the solutions to the p — A y transport equations are, as it is to be expected, polynomial in 
T. In complete analogy to the p — Z case one has that some extra conditions need to be satisfied 
in order to be able to solve the p = 4 Bianchi transport equations. More precisely, one has the 
following 

Theorem 8. Let assumptions 0-5 together with condition \84\j . Then thep — 4 Bianchi constraint 
equations are satisfied if and only if 

-^4,-l;6,3 = — io,-l;6,3- (85) 

Furthermore the functions aj^4-2q.qiT) , q = j, . . . ,A q ^ 2 are polynomial in t, while 

a,,4;4,2(r) = T4;4,2 [{I - rf-^ ln(l - t)Vj (t) + (1 + Tf+' ln(l + T)Vi-j (r) j + Q(r), (86) 
where 

T4;4,2 = ^3;4,2 ~ ^3;4,2- (87) 
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Remark 1. Again, because of the regularity conditions given in theorem if the condition 
H85|l holds then one has that 

^Te.S ~ ^3;6,3 = ^3;6,3 ~ ^3:6,3 ~ (^^) 

Remark 2. The obstruction H87(l is a result of the non-conformal flatness of our setting. It 
did not arise in the analysis carried out in |14| as the free data was chosen in such a way that 
it was automatically satisfied. It is noted that had we solved the momentum constraint in the 
way indicated in section 15.11 then the obstruction (|87|l would not have arosen — see theorem O 
As it will be seen in the sequel, analogous higher order obstruction would have, similarly, gone 
unnoticed. 



8.2.1 Solutions at order p — 5 

For the solutions of the p — 5 transport equations one has the following 

Theorem 9. Consider an initial data set satisfying the conditions 0-5. Assume, additionally, 
that the conditions i84)) , J<S,5I) also hold and the obstruction \8T)) vanishes. Then the solutions to 
the p = 5 V transport equations are polynomial in t and the transport equations implied by the 
Bianchi constraint equations at order p = 5 are satisfied if and only if 

-Zj4,0;8,4 = —Lo,o-s,4- (89) 

// the latter holds one has that aj^5-2q,qiT), j = 0, ... ,4, q — |2 — j|, . . . ,5 with q 2,3 are 
polynomial in t and 

aj,5;4,2(T) = T5;4,2((l " r)7-^"7',(T) ln(l - r) + (1 + Tf+'Vi-,{T) ln(l + r)) + Q(t), (90) 
where 

2 r. 37602 . 3099 ^ . 

Jt 5;4,2 = 18W2;4,2"l + "l^2;4,2 + mj V 677r-Lo,-2;4,2 

1023 , \ 224 



199 

Furthermore, 



\/6m(^Lo^-i;4,2 - -^4,-1:4,2) + j^V6(^Lo, 0-4^2 - -^4,0:4,2) • (91) 



a,,5;6.3(T) = T5;6,3 (^(1 " t)^"^' ln(l - t)7',+i + (1 + Tf+^ ln(l + t)7'5-,(t)) + Q(t), 



(92) 



where 

Tf"5:6,3 = ^4;6,3 ~ ^^24:6,3- (93) 

Remark 1. The expression (I91II is the natural generalisation of the quadrupolar obstructions 
obtained in paper I and in |15[ I14j. In particular, in order to make the comparison with the 
results of more transparent we note that the calculations leading to lemma 4 in section 3.2 
show that 

r-2;4,2 = 16^6^3+^2, (94) 

where 

r = r2;4,2r4V + 0(p'), 
and r is the Ricci scalar of the conformal metric as discussed in section 3. Thus, in the time 
symmetric case — i.e. if J = io.-2:4,2 = io.-i:4,2 ~ i4,-i;4.2 — -^o,0;4.2 — ^4.0:4, 2 = — the 
quadrupolar obstruction (|91|l reduces to 

^5:4, 2 = 18W2;4,2»7^^ + ™^^2:4,2- (95) 

Remark 2. Note that the first three terms in H91() are manifestly real, while the remaining 
three are — due to the reality conditions and the condition (|84fl in theorem 7 — pure imaginary. 
Thus, the vanishing of the quantity T5:4^2 implies, on the one hand, that 



8V6,, , 2089 . 

m«72;4,2 = 9-^3:4.2 - -JggJ : (96) 
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and on the other hand 



1023m(Lo,_i;4,2 - i4,-i;4,2) + 224(Lo,o;4,2 - i4,0;4.2) " 3099m^ Lo,-2-4.2 = 0. (97) 

Remark 3. As m the case of the quadrupolar obstruction (p?7jl . the octupolar obstruction 
is purely due to the non-conformal flatness of the set up and did not arise in the analysis given in 
[T^ as the initial data considered in that case was constructed so as to satisfy it automatically. 

8.3 Solutions at order p = 6 

The expansions at order p = 6 are both quantitatively and qualitatively more complex than those 
of the previous orders. Firstly, we discuss which conditions are needed to ensure the existence of 
solutions. These are given in the next 

Lemma 6. Consider initial data satisfying the assumptions 0-5, and such that that the conditions 
\8Ji\j , \85}) and hold. Furthermore, assume that the initial data is such that the obstructions 
\81\j , \91\) and \9'^ vanish. Then the solutions to the p — 6 v transport equations are polynomial 
in T and the p = 6 transport equations implied by the Bianchi constraint equations are satisfied if 
and only if 

io,l;10,5 = ^-^4,1;10,5- (98) 

The structure of the solutions to the p ^ 6 Bianchi transport equations is described in the 
following theorem. 

Theorem 10. Consider initial data sets satisfying the assumptions 0-5, and such that the con- 
ditions i84)) , 1^8 5\) . 1^8 9\) and 198\) hold. Furthermore, assume that the initial data is such that 
the obstructions i8T)) , 191\) and 1^9'^) vanish. Then the coefficients aj.6;2ij,ij(''"), with j = 0, . . . ,4, 
q = |2 — j'l, . . . , 6, q =/= 2, q =/= 3 have polynomial dependence in r. On the other hand, 

a,,6;4,2 = T+,.,{1 - Tf-'Vj{T) ln(l - r) + T-^.^^^^{1 + rf+'V^-.ir) ln(l + r) + Q(r), (99) 

where the obstructions are given by 
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58 77^,/, r \ /7T / 3282401 3648769 ^ 
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\/6^L4^1:4^2 + io,l;4,2^ 



144 

u 

7722 . „ 2198208 , . 20817 



T^6;4,2 = ^^mJQs + gQg^"" mV^ - -^AmJ^ 

62691 r- . / ^ ^ \ n: 0/7559711 8613019 



VbAm[ -L4, -1:4.2 - -t^O, -1:4,2 ^ v6to —TTTTTTTT ^0,-1:4.2 + -rTTTTTTTTr -'^4,-1:4 



r-J^ r \ /7T / 3282401 3648769 ^ 

V oA La 0-4 2 ~ n-4 2 + v om Ln n-4 2 H Fa n- 

V u.u.4,2y -TV ^ 27090 "'"'^'^ ^ 27090 ' ' 

V& \ Fa 1-4^2 + ^0. 1:4,2 



4816 V / V 126420 ' ' ' 126420 

58 
+ 43 
144 

11500 r- . , 2985 ^ , 18 , , , , 

- -^V6Amr3+ ,2 - ^V6™r4+ ^2 + ^^6(^5^4,2 - ,2)- (loob) 

Moreover, 

a,,6;6,3(T) = T6;6,3((l - rfV.+iir) ln(l - r) + (1 + rf+'V^-^ir) ln(l + r)) + Q(r), (101) 
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where 



^, . 565753248 

JL 6:6,3 = 24to u;3.6,3 ^777^7 «™ J 

82585 

36399 r 2 ( \ 1904 ^ / 

H 330"^^*"^ -^4,-1,4,2 - io,-l;4,2 I + V 6im J I L4,0;4,2 - -^-0,0;4,2 

--^^V30mLo.-i,6,3 + "^^^"^^ ^-^0,0;6,3 " -^4,0;6,3^ + %/30 ^Loa;6,3 " i4,l;6,3 

233216 r ^ ^ A- , . 
^V6*mjy3+ ^2 + 6\/30y+6,3. (102) 

Remark 1. The most noticeable feature of the order p = 6 quadrupolar obstructions (|l()()a(l 
and (|100b|l is that they are time asymmetric. That is, the vanishing of one of them does not 
imply the vanishing of the other at least if J 7^ — this feature had aheady appeared in the 
conformally flat case discussed in paper I. However, it is noted that at least to this order the 
conformal flatness plays no role in the time asymmetric behaviour. 

Remark 2. If both ^X-^^ 2 ^'^^ """6-4 2 vanish then the expressions HlOOafl and (|100b|l can be 
used to solve for, say, io,-i;4,2 and i4,_i;4,2 in a similar way to how it was done in paper 1. 



9 Conjecture about the general structure of the solutions 

The results of the CA calculations carried out in the previous section exhibit a well defined 
structure. From here and together with the results of part I, it is not hard to conjecture what 
the behaviour of the solutions for a generic order p should be. 

Firstly, we note that the conditions 184|) . (|85|) . (|89|) and H98|) arise, arguably, as a side effect of 
the way we have calculated expansions consistent with the momentum constraint Ij9b|l and that 
strictly speaking we have not found solutions. It is quite likely that the use of explicit solutions to 
the momentum constraint — as it is possible in the axially symmetric case, see section|5] — would 
get rid of such conditions. Therefore, conditions of this sort will be disregarded in the subsequent 
discussion. Moreover, following the general trend of the paper, we limit ourselves to the axially 
symmetric case. The generalisation of the discussion given in the sequel for non-axially symmetric 
initial data sets should be direct. 

As with the expansions carried out in section 5, we restrict ourselves to a class of initial data 
sets satisfying assumptions 0-5. Now, assuming that for a given integer > 7 all lower order 
obstructions vanish — i.e. those arising from the p = 4, . . . — 1 Bianchi constraint equations — 
it is conjectured that the solutions to the Bianchi transport equations are such that: 

(i) The coefficients 

aj,p;2p,p(r), aj,p;2(p-i),p-i(r), j = 0, ...,4, (103) 

associated with the harmonics T,2.p^^j^^ and T^f^^_^-^ respectively, are polynomials 

in T. 

(ii) The coefficients aj^p-2{p-2).p-2 accompanying the harmonics T'2(p-2)'' ^p-2+j ^^"^ form 

aj,p;2{p-2),p-2{T) = Tp.2(p_2),p-2(^(1 - r)P+^"-' ln(l - T)Vp-4.+j{T) 

+ (1 + Tf~^~^ ln(l + T)Pp_, (r)^ + Q(t), (104) 
where the obstruction Tp.2(p_2).p-2 is given by 

Trp:2(p-2),p-2 = ■2'p_i;2(p-2),p-2 ~ "^p- l;2(p-2) ,p-2 ' (^^S) 
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(iii) The coefficients aj-,p;2(p-3),p-3 which go together with the harmonics T'2(p-3)'' 
the form 

aj,p;2{p-3),p-3{T) = Tp.2(p_3),p-3(^(1 - t)p+^-^ ln(l - T)-Pp_5+j(T") 

+ (1 + t)p-^+' ln(l + r)Pp-i-, (r) j + Q(r), (106) 

where Tp.2(p_3),p_3 is a time symmetric obstruction similar in structure to those given by 
(|91|l and pU2|l . and should be a quantity of 2^^'^-polar nature. That is, it should be possible 
to write it as a linear combination of quantities quadrupolar nature like 

''«p-3;2(p-3),p-3, ^p-2:2(p-3),p-3' (^4,s;2(p-3),p-3 " io,s;2(p-3) ,p-3) , 

and 2^~^-polar products of lower multipole quantities such as 

^, '^■^p-3;2(p-4),p-4' '^(-^4,s;2(p-4),p-4 - io,s;2(p-4),p-4)- 

(iv) The coefficients aj^p^2q.q for g = 2, . . . ,p — 4 have also logarithmic dependence, being of the 
form 

a,,p;2,,,(r) = T+2,,,(l - r)f+2-^- ln(l - r)P,+,_2(r) 

+^;;2qA^ + rr-'^'Vq+2-,{r) (107) 

where 'Yp-2q q ^-nd '^p-2q q obstructions of 2'-polar nature. Furthermore, 

'^r,2q,q = ^ ^p;2<;,<; 0; 

that is, the obstructions are time asymmetric. 

(v) The coefficients associated with harmonics of the form and Tq*q, that is aj^p;2,i(T), 
j = 1,2,3 and a2,p;o,o(''') are polynomials in r. 



Providing a proof of the above structure would constitute a remarkable feat. As things stand 
now, it would require a much deeper understanding of the group theoretical structures underlying 
the cylinder at spatial infinity, X. 



10 Conclusions 

The discussion presented in the current article is the natural extension of the analyses carried 
out in jl5[ 1141 [T7] . The objective at the start of this programme was to provide evidence for a 
conjecture presented by H. Friedrich in 9 . As discussed at length elsewhere, this condition stated 
that only mild asymptotic conditions on the initial data were required in order to guarantee the 
existence of developments with a smooth null infinity. The conjectured asymptotic requirements 
were essentially that condition (i) in|Hlholds to all orders. As shown in 15 this conjecture resulted 
too optimistic. The results in EI have suggested that Penrose's notion of asymptotic 

simplicity imposes stringent restrictions on the behaviour of the gravitational field in the so- 
called region of spacetime "close to null and spatial infinity" . Namely, that the field should be 
static/stationary in this region — a startling rigidity result if proved to be true. 

The basis of this new conjecture was the discovery of a hierarchy of so-called obstructions 
to the smoothness of null infinity. Roughly, these obstructions are quantities which can be fully 
determined in terms of initial data sets and whose vanishing eliminates some particular singular 
logarithmic behaviour at the sets where null infinity "touches" spatial infinity. Some of these 
obstructions can be fully written in terms of the freely specifiable data — like the ones appearing 
in section (ii) and (iv) of section 9 — while others involve the coefficients Wp-2q,q which arise 
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from the expansions of the conformal factor of the initial hypersurface "d. The latter class of 
obstructions is much more complex, as it is not clear — at least at first sight — how can one 
construct data for which these obstructions vanish; the reason being that the coefficients Wp-2q,q 
are not freely specifiable. They are part of the solution to an elliptic equation (at least if one uses 
the conformal method to solve the constraints), and contain information of global nature. They 
are, in principle, (complicated) functions of the free data. 

In paper I, it was shown that the assumption regarding obstructions up to a certain order 
(p = 8) vanish leads naturally to initial data sets which are asymptotically Schwarzschildean. 
Ideally, in the present paper, one would like to repeat the same type of procedure with the 
obstructions we have been calculated, however they seem to be insufficient to be able to extract 
some type of statement concerning asymptotic stationarity. The way stationary solutions fit in 
the whole scheme, will be the subject of further study. 

Finally, we would like to make a remark concerning a point that has been implicit through the 
whole programme carried out in |15[I14[[T7| and being concluded in the present article. Namely, 
that the hierarchy of obstructions to the smoothness we have been discussed are also obstructions 
to the peeling behaviour of the spacetimes, at least in what it concerns the lower order obstructions: 
Tf"4;4;2j '^5-4,2, T5;6,3- Moreover, there are reasons to believe that the constancy of the Newman- 
Penrose constants would involve the obstructions arising at order p = 6 — i.e. '^6;6,3 and 

T6;8,4- The first steps to disentangle this connection have been sketched in ^^l- A more complete 
discussion of this connection will be given elsewhere. 
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